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The contribution to the free energy for a film of liquid of thickness h on a solid surface, due to 
the interactions between the solid-liquid and liquid-gas interfaces is given by the binding potential, 
g(h). The precise form of g(h) determines whether or not the liquid wets the surface. Note that 
differentiating g(h) gives the Derjaguin or disjoining pressure. We develop a microscopic density 
functional theory (DFT) based method for calculating g(h), allowing us to relate the form of g(h) 
to the nature of the molecular interactions in the system. We present results based on using a 
simple lattice gas model, to demonstrate the procedure. In order to describe the static and dynamic 
behaviour of non-uniform liquid films and drops on surfaces, a mesoscopic free energy based on g(h) 
is often used. We calculate such equilibrium film height profiles and also directly calculate using 
DFT the corresponding density profiles for liquid drops on surfaces. Comparing quantities such as 
the contact angle and also the shape of the drops, we find good agreement between the two methods. 

We also study in detail the effect on g(h) of truncating the range of the dispersion forces, both those 
between the fluid molecules and those between the fluid and wall. We find that truncating can have 
a significant effect on g(h) and the associated wetting behaviour of the fluid. 


I. INTRODUCTION 


by minimising the free energy functional 


The wetting of a substrate by a fluid is an impor¬ 
tant physical process and understanding such behaviour 
is crucial in a variety of fields from industrial processes 
such as lubrication and painting to biological applica¬ 
tions such as tear films in the eyes or m ucus linings in the 
lungs. The wetting behaviour of a fluid 1 ^ is determined 
by the manner in which the atoms or molecules within 
the fluid interact with each other and with those forming 
the substrate. Determining the macroscopic fluid prop¬ 
erties, wetting behaviour and thermodynamics, starting 
from an understanding of the (microscopic) molecular in¬ 
teractions is one of the cornerstone problems in liquid 
state science. 2 

On the macroscopic scale, the wetting behaviour of 
a fluid in contact with a solid substrate can be charac¬ 
terised by the contact angle that a liquid drop makes with 
that substrate. Three regimes of wetting behaviour can 
be identified: complete wetting, partial wetting and non 
wetting, these three states are defined by contact angles 
of 6 = 0°, 0° < 9 < 180° and 9 = 180° respectively.^ 
Surface tension forces arise from interfaces in the fluid 
and these can be related to the contact angle by Young’s 
equation^ 

COS»= ---, (1) 

7 lg 

where 7 i g , 7 w i and 7 wg are the liquid-gas, wall-liquid and 
wall-gas surface tensions respectively. 

The effective interface Hamiltonian (IH) model} 3 also 
referred to as the interface free energy model, describes 
the height profile of a mesoscopic liquid film on a sub¬ 
strate. One can find the equilibrium shape of a droplet 


F[h}= J [<?(/>)+W1 + W 


dx, 


( 2 ) 


where h(x) is the liquid film thickness at some point x 
on the substrate and g(h) is the binding potential, which 
is also referred to as the effective interface potentialP® 
g(h) is a restricted free energy, i.e. the free energy sub¬ 
ject to the constraint that the thickness of the liquid layer 
adsorbed on the surfaces is h. A good discussion on the 
subject of restricted free energies can be found in Ref. j5] 
The binding potential describes the interaction between 
two interfaces and is related to the disjoining pressure 
n = —dg/dh. In the IH model, the binding potential is 
often approximated by simple expressions that give the 
qualitatively correct behaviour. A common example of 
such an approximation would be an asymptotic expan¬ 
sion which is valid only for larger film thicknesses (cf. 
Sec. 0 - This paper sets out a method to directly cal¬ 
culate the binding potential from a microsc opic basis, 
namely via density functional theory (DFT) rb 0 ' ' 15 As 
input, the method takes the interactions between parti¬ 
cles in the fluid and also the forces on the fluid particles 
due to any external fields, for example that due to the 
wall of a container or a surface on which the liquid is 
deposited. This calculation yields an expression for the 
binding potential that is valid for all film thicknesses. 
The term “particle” is used here generically to refer to 
the atoms/molecules/colloids in the fluid, depending on 
the exact system under study. Note that it is also pos¬ 
sible to calculate g(h) from computer simulations - see 

Refs.QMSI 

The IH model is particularly useful due to its appli¬ 
cation in dynamical studies. By making the assumption 
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FIG. 1: Two different descriptions of a liquid drop on a sur¬ 
face: (a) a height profile calculated via Eq. 0 from a meso¬ 
scopic free energy (cf. Eq 0 ) and (b) a density profile, which 
gives the fluid number density at a distance 2 above the sur¬ 
face. These are both for a fluid with /3e = 0.9 and f3e w = 0.6 
(see Sec. Ill for further details). 


of small surface gradients and contact angles, Eq. © re¬ 
duces to 


m 


#) + f(vft) 2 


dx 


(3) 


where we have omitted a constant contribution. Equa¬ 
tion 0 can then be emp loyed in the thin him (or long¬ 
wave) evolution equation^ 25 which describes the time 
evolution of a thin film of liquid on a fl at sol id substrate. 
In gradient dynamics form it is writterP§l2Z 



Q(h)V 


SF[h]' 

Sh 


(4) 


where Q{h) is a mobility factor that depends on the film 
thickness h(x, t). Equation @ may be derived by mak¬ 
ing a long wave appr oximation in the governing Navier- 
Stokes equation ] 25 * 28 ^ There are many applications of this 
equation to model different situations. Steady state solu¬ 
tions, where dh/dt = 0, such as drop profiles, are found 
by minimising the free energy, Eq. 0 , with respect to 
the film height profile, subject to a volume constraint. 
More specifically, it amounts to solving 



(5) 


where a is a Lagrange multiplier stemming from the con¬ 
straint 


J dx /i(x) = Vo, 


( 6 ) 


where Vq is a specified drop volume. A typical drop 
profile resulting from such a calculation can be seen in 
Fig0a). Note the very thin non-zero height ‘precursor’ 
film that is present to the left and right of the drop. 

A fully microscopic description has statistical mechan¬ 
ics as its basis. Statistical mechanics calculates an aver¬ 
age over an ensemble of all possible states of the system, 


i.e. it averages over all possible configurations of the par¬ 
ticles. This average leads to determining the fluid one 
body density profile p( r) , which represents the likelihood 
of finding a particle at a given point r in the systemP 
This statistical mechanical point of view is the basis for 
DFT. From DFT, the grand potential, f2, of the system 
is calculated and the equilibrium density profile, p( r), 
which minimises Cl, can be found. A typical equilibrium 
density profile of a liquid drop on a solid substrate is 
displayed in Fig. 0b). Other examples of drop profiles 
calculated using DFT can be found in Refs. [2M551 Note 
that by identifying the surface of the liquid as the surface 
where the density equals a specified value, p mt , where 
Pi > Pint > Pg and where pi and p g are the coexist¬ 
ing liquid and gas densities of the fluid, the description 
can be further reduced to obtain a film height profile 
very similar to that displayed in Fig. 0a). One possible 
choice is to choose pi nt = ( pi + p g )/2 . Alternatively, by 
integrating in the 2 direction over such a density profile 
we can obtain a drop height profile. Here, we define the 
height of the liquid film on a substrate as the adsorption 
divided by the liquid-gas bulk density difference 

^ (| = £Ww)Al£. (7) 

Pi - Pg 


It is worth noting that using DFT it is just as easy to 
compute the profile for a liquid droplet that makes a 
contact angle with the substrate that is greater than 90°, 
than one with 6 < 90° . It is significantly more difficult to 
find drop profiles for 9 > 90° in the mesoscopic approach 
because for these contact angles, Eqs. (0 and 0 cannot 
be employed. 

The remainder of this paper is set out as follows: in 
Sec. |0 the binding potential and the procedure for calcu¬ 
lating it are discussed. A simple DFT, for the lattice-gas 
model, is presented in Sec.[nl]that is used to demonstrate 
the procedure. The dependence of the fluid behaviour 
on the particle interactions is discussed in Sec. IV In 
particular, it is shown that truncating the range of the 
dispersion interactions between the particles has a pro¬ 
found effect on the binding potential and interfacial phase 
behaviour. The method of fitting a function to the cal¬ 
culated data is given in Sec. [Vj followed by the results of 
passing the binding potential from the lattice-gas model 
to the thin film IH model in Sec. Eg Finally, conclusions 
are drawn in Sec. IVII1 


II. THE BINDING POTENTIAL 

For any fluid, coexistence between the liquid and gas 
phases occurs when the temperature T, chemical poten¬ 
tial p, and pressure p of the two phases are equal. It then 
follows, for a given volume V, that the bulk grand free 
energy, Cl = —pV, of either phase occupying the same 
volume is also equal. For a system where the liquid and 
gas phases of a fluid exist together at the point of liquid- 
gas coexistence, any excess, over bulk, contributions to 
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the free energy of the system must stem solely from the 
interface that forms between the two phases. This ex¬ 
cess grand potential per unit area of the interface defines 
the surface tension between those two phases. In this 
case, it is the liquid-gas surface tension 7 i g . Now con¬ 
sider a system with chemical potential p = p CO ex, the 
value at coexistence, where a film of liquid separates the 
bulk gas from a solid surface (cf. Fig. [ 2 ]). The excess free 
energy now consists of the sum of the two interfacial ten¬ 
sions and also the interaction between the two interfaces. 
The excess grand potential in such a system is given by 
f l ex (h) = f] + pV = u ex (h)A, where A is the area of the 
interface and 


w e x (h) = 7 w i + 7 i g + g(h). ( 8 ) 


h is the liquid film thickness, 7 ^ and 7 i g are the wall- 
liquid and liquid-gas surface tensions respectively.^ The 
final term g(h) is the binding potential, which gives the 
contribution to the free energy from the interaction be¬ 
tween the two interfaces. This has the property that as 
h —¥ 00 , g(h) —>• 0. The absolute minimum of the grand 
potential defines the equilibrium film thickness h. 

In the case of liquid droplets surrounded by their 
vapour on a solid substrate, the absolute minimum of 
the binding potential is directly related to the equilib¬ 
rium contact angle of the droppSEH 


6 = cos 1 



(9) 


where <?(/i 0 ) is the value at the minimum of the binding 
potential. This corresponds directly to Young’s equation, 
Eq. ([l]); the absolute minimum of the binding potential 
corresponds to the equilibrium state of the system and 
gives the equilibrium excess grand potential, i.e. the wall- 
gas surface tension: 


uj(h 0 ) = j w i + 7 ; g + g(h 0 ) = ^ wg . (10) 

Replacing 7 wg in Eq. 0 with this expression leads to 
Eq. ((£]). 

Equation ([ 8 ]) is given above as a function of the film 
height h. From a microscopic viewpoint it is often more 
convenient to use the adsorption of the fluid as the or¬ 
der parameter characterising the fluid at the interface, 
instead of the film thickness. The total adsorption is 
readily calculated from the fluid density profile as 

r = j J dr(p(r) - p b ). (11) 

pb is the bulk fluid density, which in the cases considered 
here is the gas density p g . We may also define a local 
adsorption [c.f. Eq. 0] as follows: 

/»00 

F {x,y)= l [p(x, y, z) — p g \ dz, (12) 

Jo 

so that when the definition in Eq. 0 for the film height 
h is used, we have T = h(pi — p g ). For other definitions 


gas 




liquid 

l 

///////////////////////// 

solid 


FIG. 2: A schematic of the system: a film of liquid of thickness 
h separating a semi-infinite volume of gas from a solid surface. 


of the film height h, then 


r Khfa-Pg). (13) 

The dependence on how precisely h is defined becomes 
negligible, for large h. Thus, Eq. Q and also the binding 
potential g may be given as a function of the adsorption. 
Note, however, that negative adsorptions are possible, 
e.g. at a purely repulsive wall, and in such a situation 
describing the liquid film at the wall via a film thickness 
h becomes meaningless since the quantity h defined in 
Eq. ([?) then becomes negative. 

For a system with chemical potential p = p CO ex + 
i.e. off-coexistence, then there is an additional contribu¬ 
tion T5p, that must be added to the right hand side of 
Eq. Together with Eq. 0§, this gives 

u ex {h) « 7 w i + 7 i g + (pi - p g )h5p + g(h). (14) 


Often, only asymptotic forms of binding potentials are 
used which are strictly valid only in the limit of a large 
film thicknessP There are two main asymptotic forms 
of the binding potential that are considered, the choice 
of which depends on the assumed particle interactions 
and the range of those interactions. For van der Waals’ 
(dispersion) i ntera ctions, the following asymptotic form 

is appropriateP2§HS3 

9(h) = p + ^3 H , (15) 

The equivalent disjoining pressure i d 24 * 25 l 

n W = P + £ + -- (16) 


Under certain approximations the asymptotic behaviour 
shown in Eq. (151 can be calculated analytically includ¬ 
ing the values of the coefficients a and 6 , and how they 
depend on the temperature. 

With only short ranged interactions between particles, 
the binding potential can be expressed asymptotically 

a j3l7l8l41l42l 


g(h) = aexp(—/i/£) + &exp(—2/i/£) + • • • . (17) 
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FIG. 3: A series of binding potentials for a fluid with fixed 
inverse temperature /3e = 0.9 against a solid substrate with 
varying attraction strength /3e w . A change in wetting be¬ 
haviour, from non-wetting to wetting occurs, as indicated by 
the change in the position of the minimum in g(T) from a 
low finite adsorption to a large (T —> oo) adsorption, as /3e w 
increases. 


The length £ is the bulk correlation length in the liquid 
phase at the interface. These asymptotic forms, trun¬ 
cated after a few terms, are used frequently throughout 
the literature (independently or as combinations of the 
two form!] 1 1 25 I 4 S1M1) even though they are only strictly 
valid for thick liquid films and can not describe the bind¬ 
ing potential as h —> 0. To describe the small h be¬ 
haviour, the full form of g{h) is required. A fully mi¬ 
croscopic theory, that describes the fluid structure at the 
wall, is required to obtain this. 


Fig-! shows binding potentials for various values of 
a parameter e w , which determines whether or not the 
fluid wets the wall. These binding potentials are calcu¬ 
lated using the microscopic DFT-based approach that is 
introduced below (cf. Sec. III). The parameter e w char¬ 
acterises the interaction strength between the wall and 
the fluid particles. The full form of the long range (al¬ 
gebraic) dispersion interactions between the particles are 
included and so for large T, g(T) has the asymptotic form 
given in Eq. (151. In Fig. [3] we see that for small values of 
e w (weakly attracting, solvophobic wall) the global min¬ 
imum in g is at a small value of T, i.e. the liquid does 
not wet the wall. As e w is increased, there is a first order 
wetting transition when /3e w « 0.74 and for j3e w > 0.74 
the global minimum in g(r) is at F —» oo, i.e. there is a 
macroscopically thick film of liquid on the wall at coex¬ 
istence g = g-coex- 


DFT gives a route by which the free energy may be 
calculated, taking into account the microscopic structure 
of the fluid at the wall and the interactions of the particles 
within it. The equilibrium state of the system is found 


by minimising the grand potential functional 121MEI 

^[P(r)] = F[p(?)\ + J drV ext {r)p{r)-g J drp(r), (18) 


where T is the intrinsic Helmholtz free energy functional, 
g is the chemical potential and V ex t is the external poten¬ 
tial. The equilibrium density profile, p(r), is that which 
satisfies the Euler-Lagrange equation 


(5f2[p] 

Sp(r) 


= 0 . 


(19) 


Once an approximation for T is specified, this equation is 
solved numerically via a scheme in which an initial den¬ 
sity profile is supplied and then iterat ed until a specified 
convergence criterion is reached] 45 ! 46 ! 

In order to calculate the binding potential as a function 
of the adsorption using DFT, it is required to evaluate 
the free energy for any specified F, i.e. in addition to 
the equilibrium profile we require other non-equilibrium 
profiles for a range of values of the adsorption T. By us¬ 
ing the procedure developed and justified in Ref. [47] (see 
also Ref. [48|, the excess density (p(r) — pb) is normalised 
at each iterative step, which is equivalent to including 
a (fictitious) additional effective external potential that 
stabilises a wetting film of the desired thickness with the 
specified value of T. Using this normalisation procedure, 
the excess grand potential can be obtained for a range 
of values of T. The binding potential is then given via 
Eq. ([8]) by subtracting the values of the solid-liquid and 
liquid-gas surface tensions. 


III. A MICROSCOPIC MODEL 

The method outlined above for calculating the bind¬ 
ing potential is valid for any DFT model. To illustrate 
the procedure, a simple approximate DFT for a lattice- 
gas, is used. The model is only briefly described here: a 
full description and derivation can be found in Ref. [45] 
The system is discretised by a cubic lattice and the fluid 
particles, assumed to all be identical and spherical, oc¬ 
cupy only one cell each on the lattice. The diameter of 
each particle, and the width of each cell, is a and there 
are M = M t MjM k lattice cells and N < M fluid par¬ 
ticles in the system. A point in the lattice is denoted 
i = [i,j, k ) and Mi, Mj, and Mj. are the number of cells 
in the i, j and k directions, respectively. Note i, j and k 
are integers. Any configuration of particles in the system 
can then be described by the set of occupation numbers, 
{rii}, where rq = 1 if there is a particle in cell i and 
Hi = 0 otherwise. The average of an ensemble of all such 
systems may be taken and then the average occupation 
number of each cell is described by the density pi = (rii). 
It may be assumed that the equilibrium density profile 
only varies in one direction, or, as for the calculations 
below, it is assumed that the density profile is invariant 
in the third dimension, so only a two-dimensional (2D) 
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slice of the system needs to be studied. Here, the z-axis 
direction, indexed by k, is assumed to be invariant. The 
energy of the full three-dimensional system is given by 
the Hamiltonian 

M M 

E = - oEE^j +Y^ nh ( 2 °) 

i=l j/i i=l 


where V\ is the external potential and eyj is a Lennard- 
Jones-like pair potential between particles at lattice sites 
i and j: 


£i.j = v( rii ) 


-e/nj 6 for nj > a, 
oo for n j < a, 


( 21 ) 


where nj is the distance between a pair of particles lo¬ 
cated at lattice sites i and j. If one assumes that the sys¬ 
tem is invariant along the direction indexed by k, then 
a mean-field approximation for the internal energy U is 
the following 2D sum: 


A?2 d 


M 2 d. 


U — (E) — FI e io,joPioPjo + FI IdoPio; (22) 


io=l jo^io 


io = l 


where sums over the M 2 d = MiMj sites in the 

2D lattice plane where k = 0 and io = (i,j, 0). The third 
(invariant) dimension is now accounted for in the interac¬ 
tion weights and V\. The sums are written explicitly 
in Eqs. (B1) to (B5) in Appendix [B| these equations also 
give the weights when the range of particle interactions 
is truncated. The grand potential of this system can be 
approximated as follows®! 


n({pi 0 }) = U-TS-nN, 

M 2 d 

= k B T Y [Pi 0 hi(Pio) + (1 - Pi 0 ) M 1 - Pi 0 )] 

io=l 

2 M2d M2d 

~ 9 5 Z F/ e io,jo/ 3 >oPjo + F/ ^‘ 0(^0 — p)i 
" io=ljo^io io = l 

(23) 

where k B is Boltzmann’s constant and S is the entropy 65 . 
Included implicitly in the above equation is the particle 
diameter a = 1 . 

The set of lattice densities that describes the system 
at equilibrium is found by solving [c.f. Eq. (19)]: 


d Pi ’ 


(24) 


for every lattice site i. From Eqs. (23) and (24) we obtain 

. MiMj 

k B T In f - j - Y e i.jPj + Vi ~ P = 0. (25) 

V pi J j=1 


Solving this coupled set of equations (25) gives the equi¬ 
librium fluid profile {pi}. However, in order to obtain 


the non-equilibrium profile with specified adsorption T^, 
the iterative method developed in Ref.[47]must be used. 
This is equivalent to solving the set of coupled equations 
[c.f. Eq. 


k B T In 


Pi 


1- Pi 


Mi Mi 


-- Y ^jPj + Vi + Vf-p = (l (26) 


j=i 


where V] 6 ® is the fictitious additional potential mentioned 
at the end of Sec. [h] It has the property that E ; eff -» 0 
as i —¥ oo, far from the wall. This additional potential 
stabilises a film of liquid with the specified adsorption 
against the wall. Note that V[ eff is a-priori unknown, 
it is calculated on-the-fly self consistently as part of the 
minimisation algorithm. 47 ! This is done by re-normalising 
the density profile during every iteration of the algorithm 
by replacing the value of the density with 

pr w = ^ L (pf d -p b )+p b , ( 27 ) 


where T^ is the desired adsorption and Tq^ is the ad¬ 
sorption corresponding to the density profile {p? ld }, ob¬ 
tained from iterating Eq. (25). More details about this 
algorithm and its properties can be found in Ref. [47] 

Note that this procedure does not require the bulk 
phase to be at coexistence, although in all the results 
presented here it is at coexistence. As the bulk fluid 
state point (p,T) is varied, the form of the restricted 
free energy p(T) also changes as well as the bulk den¬ 
sities and, in consequence, the interface tensions. Note 
also that one can vary the adsorption T by the standard 
method of varying the value of the chemical potential 
p. Calculations (not presented here) show that the main 
difference between results from our method for varying 
the film thickness via a fictitious external potential, with 
results from varying it by changing p are to be seen when 
the adsorption is small. This is particularly so when the 
fluid wets the wall, because in this case, to obtain a small 
adsorbed film height by varying the chemical potential 
requires a large shift in the value away from that at co¬ 
existence. 

The bulk fluid phase diagram is displayed in Fig. [4] 
For details about how this is calculated see e.g. Refs. 
I45|49l The binodal and spinodal are both displayed. The 
binodal gives the densities of the coexisting gas and liquid 
states. Within the spinodal curve, the uniform fluid is 
unstable and spontaneous demixing occurs. The bulk 
critical point is at pa 3 = 0.5 and temperature k B T/e = 
1.5. 


DFT is a statistical mechanical theory - i.e. in prin¬ 
ciple it should give the ensemble average density profile 
of the fluid. A statistical description of a fluid confined 
in an external potential should yield an (ensemble av¬ 
erage) density profile with the same symmetry as that 
potential) 10 ! 5 ^ ! Thus, for a planar wall, the equilibrium 
density profiles only vary with the distance from the 
wall. Fig. [5] shows typical examples of such profiles for 
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FIG. 4: The phase diagram for the lattice fluid in the temper¬ 
ature kgT/e versus density per 3 plane. The solid (red) line is 
the binodal curve and the dashed (blue) line is the spinodal 
curve. 


the lattice-gas model and of the corresponding points on 
the binding potential that they represent. This series of 
density profiles range from small (including negative) val¬ 
ues of the adsorption to large values, where the profiles 
indicate that there is a thick film of liquid at the wall. 

The corresponding fictitious potentials are also dis¬ 
played in Fig.[5j At the minimum in <?(r), the fictitious 
potential VJ eff is, of course, zero, because this corresponds 
to the equilibrium state. Moving away to either side of 
this minimum, V( eff increases rapidly and the largest mag¬ 
nitude potentials are observed for low adsorptions T, for 
values of T where the gradient in g(T) is largest. In Fig. [5] 
the density profile (and associated potential Vj efI ) corre¬ 
sponding to Ter 2 = 0.304 is very close to at equilibrium 
and so Id eff is very small everywhere. In contrast, the 
fictitious potential for Fa 2 = 0.004 is much larger, as the 
gradient of g(T) at this point is also large. For larger ad¬ 
sorption values, in the tail of the binding potential, Id 6 ® 
can be either weakly attractive or weakly repulsive. This 
is due to the oscillations in p(T) stemming from the fact 
that we are dealing with a lattice model. In a continuum 
DFT model these oscillations decrease in amplitude as T 
increases or are entirely absent, depending on the fluid 
state point. In the subsection below, we discuss this issue 
further. 

Note that even when the typical microstates of the 
system consist of liquid drops on the surface, after per¬ 
forming a statistical average over all states, the resultant 
density profile should be invariant in the direction paral¬ 
lel to the substrateP3 In order to study liquid drops or 
liquid ridges (invariant in one direction along the surface) 
one must break the translational symmetry and impose 
that the centre of mass be located at a particular point or 
line on the surface. By constraining the centre of mass, 
2D drop profiles (i.e. liquid ridges) can also be found for 
the same system even though the external potential only 


varies in one direction. 

Fig.[6]displays such 2D density profiles for external wall 
potentials of varying attraction strengths. These are cal¬ 
culated by solving Eq. (25). However, the value of p is 
not at the outset imposed, instead the total adsorption 
is specified. This is done by renormalising the den¬ 
sity profile at each iteration of the algorithm via Eq. ( [27) . 
For more details about this algorithm, see Ref. 05J 

Fig. [6] shows that for weakly attracting walls (small 
e w ), it is energetically favourable for the liquid not to be 
in contact with the wall and the contact angle is large. As 
e w is increased, the contact angle decreases as the fluid 
seeks to have greater contact with the wall. As discussed 
below (see also Fig. [3]), there is a wetting transition at 
j3e w ~ 0.74 and for values of /3e w greater than this, the 
liquid spreads over the surface (9 = 0). Note that the 
lattice-gas DFT, with various different choices for e;j, 


has been used extensively to study wetting and also to 
calculate profiles for liquid drops on surfaces - see e.g. 
Refs. H7TH551 for some recent example. As far as we are 
aware, in all previous studies the adsorption on the sur¬ 
face is determined by the choice of p, either explicitly, or 
by fixing the total number of particles in the system, N. 
Since the fictitious potential that we use in our method 
to stabilise the drops is almost always rather small, the 
density profiles we obtain are actually rather similar to 
those found previously when the interaction and wall po¬ 
tentials are the same. 


A. The Binding Potential via the Sharp-Kink 
Approximation 

Some of the interfacial and wetting behaviour of the 
lattice-gas model can be elucidated in a particularly sim¬ 
ple manner by making the so-called sharp kink (SK) 
approximationPEI Specifically, the lattice gas density pro¬ 
file is taken to be 

! 0 for i < 0, 

pi for 0 < i < h, (28) 

p g for i > h, 

where h is the position of the liquid-gas interface and the 
surface of the solid wall is at the lattice site i = 0. No 
minimisation is necessary under this approximation and 
the binding potential for any given film thickness h can 
be calculated directly. For a system with truncated fluid- 
fluid interactions, the binding potential can be calculated 
as 


OO 

g{h) = (pi - p g )a 3 ^2 V, (29) 

i=h-\-1 


when h > A, the range of the fluid-fluid interactions. An 
asymptotic expansion of this sum gives (see also Ref.© : 


9 


e w n(pi - p g )a 3 



12h 3 



(30) 
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FIG. 5: In (a) we display the binding potential g( F) for a fluid with bulk gas density per 2 = 0.107 and temperature ft = 0.9 
against a planar wall with attraction strength f)t w = 0.6. In (d) we display a magnification of tail of g(T), for larger values of 
the adsorption T. Marked on p(F) are points that correspond to the density profiles displayed in (b), and the corresponding 
fictitious potentials Vf s , which are displayed in (c) and (e). The marked points have adsorption values: Ter 2 = 0.004, 0.304, 
0.6, 1.6, 3.6, 4.6, 5.6, 6.6, 7.6, 8.6, 9.6 (note that the density profiles for the final four values are not displayed in (a) and (b), 
for clarity). The fictitious potential is that which must be applied to stabilise a film of liquid with the given adsorption in an 
open (grand-canonical) system. We observe a clear relation between the gradient of p(r) and the magnitude and sign of the 
fictitious external potential. 


Truncating this series gives an approximation that is 
valid for large h. The coefficients given in Eq. (15) can 
therefore be calculated explicitly for the lattice-gas model 
under this approximation. Figure [7] shows the binding 
potential for fit = 0.9 and fie w = 0.55 on a log-log plot, 
calculated using the full lattice-gas model and comparing 
with results from the SK approximation. The analyti¬ 
cally calculated asymptotic limit in Eq. (30) is plotted to 
0 (/W 2 ), which agrees very well with a numerical evalu¬ 
ation of the full SK approximation. The numerical and 
analytic SK results begin to slightly drift apart for larger 
adsorptions. This is because in the numerical calculation 
the external potential is truncated to a range of 100 ct 
whereas the analytic calculation assumes an infinite in¬ 
teraction range. This truncation results in the numerical 
result diverging from the ~ h~ 2 decay in Eq. (301 for 
large h. 


The density profile obtained by minimising the DFT, 
in contrast to the SK approximation, incorporates a bet¬ 
ter approximation for the true shape of the liquid-gas 
interface and also the effect of this interface being close 
to the wall. Thus, for small values of T (i.e. small h), this 
approximation is far more reliable. Also, it includes the 
correct T -2 asymptotic decay for large T. Note that in 
Fig-0 the oscillations in the tail of the binding potential 
are present as a result of the system being discretised 
on a lattice. The free energy is lower when the liquid- 
gas interface is between two lattice sites, rather than on 
a lattice site. Thus, as the specified adsorption is in¬ 
creased, the liquid-gas interface moves continuously away 
from the wall which leads to oscillations in g(T). These 


oscillations he on top of the correct ~ T~ 2 asymptotic 
decay; i.e. for large h the binding potential is of the form 
g ~ e w na 5 (pi — p g )/(12h 2 ) — B cos(27r/i/cr), where the 
amplitude of the oscillations B is a small number that 
depends on the state point and manner in which the in¬ 
teractions are treated. As the range of the fluid-fluid 
interactions is increased, the amplitude of these oscilla¬ 
tions B decreases (c.f. Fig.[sj) . Thus, apart from these os¬ 
cillations, the binding potential calculated from the full 
minimisation matches up well for large T with the results 
from the SK approximation. We attempted to extract 
the coefficients for the higher terms in the expansion in 
Eq. (15) from our numerical results, since analytic ex¬ 
pressions for these existP^ However, the oscillations in¬ 
duced by the lattice in g{h) make this problematic. Note 
too that the free energy contribution from the liquid-gas 
free interface "fi g is a little different in the SK approx¬ 
imation compared to that obtained from the full min¬ 
imisation. The respective values are subtracted when 
calculating p(T) in the two methods. 


For small values of the adsorption, large differences are 
seen between the SK results for g(T) and those from the 
full minimisation. This demonstrates that the SK ap¬ 
proximation should only be used for thicker films and a 
microscopic theory for the fluid structure (DFT) is re¬ 
quired to accurately calculate a binding potential that is 
valid for thin (h < 3<r) liquid films. 
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FIG. 6: A series of 2D density profiles for a drop of the lattice- 
gas fluid with /3e = 0.9, deposited on various solid substrates 
with attraction strength parameters /3e w = 0.2, 0.5, 0.6, 0.72 
and 0.74. The fluid-fluid particle interaction range is trun¬ 
cated to L = 5. 


IV. THE INFLUENCE OF THE RANGE OF 
PARTICLE INTERACTIONS 


In Sec.[TT]the manner in which the particle interactions 
influence the (asymptotic) form of the binding poten¬ 
tials was discussed; i.e. for long range dispersion interac¬ 
tions the algebraic decay form in Eq. (15) is appropriate, 
whereas for short-range forces, the exponential form in 
Eq. (JTtJ) is the correct asymptotic form. In this section, 
the cross-over from one regime to the other is explored 
as the interaction potentials in the lattice-gas model are 
truncated at different ranges. 


When comparing systems with different interaction 
ranges, it is important that the truncation does not 
change the bulk fluid phase diagram. Otherwise, even 



FIG. 7: The binding potential for /3e = 0.9 and /3e w = 0.55 
calculated for the lattice gas model and compared with re¬ 
sults from the sharp kink approximation, displayed on a log- 
log plot. The green dashed line corresponds to the leading 
order (~ h~ 2 ) term in Eq. < |30| ) and the dotted blue line to 
a numerical evaluation of the sum in Eq. (291. The red solid 
line is the result from a constrained minimisation of the full 
functional. 


if comparing systems at the same temperature T, the 
value of (T — T c ), where T c is the bulk critical tempera¬ 
ture, would not be the same for systems with a different 
truncation range L. For the simple mean-field approxi¬ 


mation to the free energy in Eq. (23), the bulk uniform 
fluid free energy is determined by the integrated inter¬ 
action strength of the pair potential i.e. the to¬ 

tal potential that arises from the interaction of a sin¬ 
gle particle with all others within the interaction range. 
As the interaction range is adjusted, one must vary the 
value of e, the parameter governing the overall strength 
of the pair interactions, to ensure that the integrated 
interaction strength remains constant so that the bulk 
fluid phase diagram remains unchanged. All values of f3e 
quoted here are the strength of the interaction when the 
interaction range is truncated to only the nearest neigh¬ 
bour lattice sites (L = 1). 

Truncating the range over which particles in the sys¬ 
tem interact changes the overall shape of the resulting 
binding potential. In particular, if all interaction poten¬ 
tials (both fluid-fluid and wall-fluid) are truncated, then 
the tail of the binding potential decays exponentially to 
zero. If there are any long ranged (not truncated) inter¬ 
actions then, in three dimensions, the binding potential 
tail decays algebraically ~ h~ 2 . By varying the range 
of the particle interactions in the lattice-gas model, both 
of these regimes are seen and also there is a crossover 
from one to the other as the truncation range L is var¬ 
ied. These results are displayed in Figs. [8] and [9] 

Fig. © shows the binding potential for (3e = 0.8, 
/3e w = 0.5 and at different truncation ranges, L, of the 
fluid-fluid interactions. Note, however, that the inter- 
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actions between wall and fluid particles extends over the 
entire domain in all cases (the domain size is Mk = lOOer). 
As the truncation range L is increased, it becomes ener¬ 
getically favourable for the fluid not to wet the substrate 
and so a minimum in g(T) at a finite value of the ad¬ 
sorption F appears; i.e. the interfacial phase behaviour 
changes purely as a result of how the fluid-fluid particle 
interactions are modelled. Note that since the value of 
g(r) at the minimum g(To) = g(ho) decreases as L in¬ 
creases, from Eq. Q, this indicates that increasing the 
interaction range L makes the fluid less wetting and in¬ 
creases the contact angle. This shows that care should 
always be taken when modelling the interaction between 
two fluid particles: truncating the pair potentials at too 
small a distance may result in significant errors in pre¬ 
dictions for the interfacial phase behaviour. Fig. [ 8 ] also 
shows the ~ T -2 decay of g(T) as T — > oo that is present 
in all cases, due to the presence of the long ranged wall- 
fluid interactions. 

The effect of truncating the range of all interactions 
(including the wall-fluid potential) on the binding poten¬ 
tial is shown in Fig. [9] for the case when (3e = 0.8 and 
/ 3e w = 0.7. The binding potential for L = 80 is the 
longest ranged, having the slowest decay to zero as T in¬ 
creases. As L is decreased, the range of g(T) decreases. 
The inset of the figure shows the same data on a log-log 
scale which enables one to clearly see the form of the 
asymptotic decay of g(T) as T —>• oo. As expected from 
the discussion above, as the range of the interactions L 
is increased, an increasingly large portion of algebraic 
~ T ' 2 decay is present. However, it should be pointed 
out that formally speaking, it is only when L —»• oo that 
the ultimate asymptotic decay of g(T) as T — > oo changes 
from exponential to algebraic. 


V. A FITTING FUNCTION FOR THE BINDING 
POTENTIAL 


In order to take the binding potentials calculated via 
DFT in the previous section and use them with the 
mesoscale IH model, a suitable fit function is required. 
The fit function should have the same form as that ob¬ 
served in the DFT results so that it can accurately fit the 
data. Several aspects of the binding potential curves are 
particularly important for the fit function to be correct: 
At small values of the adsorption T, the binding poten¬ 
tial exhibits a minimum and a maximum when the fluid 
is non-wetting or near to the wetting transition. These 
need to be fitted well; in particular the value at the min¬ 
imum g(ho) needs to be the correct value [c.f. Eq. Q]. 
Also, the true binding potential is finite for T = 0 (un¬ 
like the approximation in Eq. (15) and other such power- 
series) and so the fit function should not diverge at T = 0 . 
Finally, when dispersion interactions are present, the fit 
function should exhibit the correct T -2 decay as T —» oo. 
Therefore, the following form for the fit function is sug- 



FIG. 8: The binding potential for the fluid with inverse 
temperature /3e = 0.8 at a wall with attraction strength 
Ptw = 0.5. The various different curves correspond to trun¬ 
cating the fluid-fluid pair interactions at different values of the 
truncation length L. The wall-fluid interactions remain fixed 
at a truncation of lOOcr. Note that we vary L in a manner 
that does not change the bulk fluid phase diagram. We see 
that the fluid is predicted to be more wetting as L becomes 
shorter; in fact, the interfacial phase behaviour changes from 
non-wetting to wetting, purely as a result of truncating the 
interaction range. The inset shows the same curves plotted on 
a log-log scale which shows that the asymptotic decay form 
for large T is the same in all cases because of the long-ranged 
wall-fluid interactions. 



r<r 2 


FIG. 9: The binding potential for the fluid with inverse 
temperature /3e = 0.9 at a wall with attraction strength 
/3e w = 0.7. For these values, the fluid wets the wall. In 
contrasts to the case in Fig. |8j here we truncate both the 
wall-fluid and fluid-fluid interactions at the same range, L. 
The different curves correspond to varying L from between 1 
to 80 particle diameters. As L is decreased, the range of (/(F) 
decreases. In the inset, the black dashed line is O.OOir -2 . As 
L is increased, the binding potential approaches this large T 
asymptotic decay form. 
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gested: 


9( r) = ^ p l-y 2 r) ]- 1 , 


(31) 


with 


P(T) = a 0 r 2 e- aii + a 2 T 2 + a 3 + a 4 r 4 + o 5 r 5 . (32) 

The rationale for this choice is as follows: For small x, 
exp(x) ~ 1 + £, and so at low adsorptions this form gives 
g(T) « A{a 0 e~ aiT + a 2 + a 3 T + ...). For high adsorp¬ 
tions this form gives g(T) ~ — AT -2 , which is the correct 
form for the asymptotic T —> oo decay. The coefficient 05 
of the highest order term must be positive. The second 
exponential within the outer exponential function [i.e. 
the exponential term with coefficient ao in P(T)] helps 
to correctly fit the minimum at low adsorptions, which is 
often asymmetric, being much steeper on the low adsorp¬ 
tion side of the minimum, compared to the other side, as 


can seen in Fig. 10 The constant a\ in Eq. (31) is usu¬ 


ally quite large and positive so that the inner exponential 
term has almost no effect on the form of g(T) on the large 
adsorption side of the minimum. 

An example of the fitting function is displayed in 
Fig. [lO] The function is plotted with the original data 
and the SK results as a comparison. The SK results do 
not describe the behaviour at small T whereas the fit 
function gives a very good approximation to the data 
over the whole range. The leading order coefficient of 
the decay of the binding potential, A in Eq. (31), can ei¬ 


ther be calculated directly using the SK approximation 
or can be found by fitting the DFT data, both methods 
give similar results. The value for A from the SK ap¬ 
proximation is used here. In Appendix A we give the fit 
parameter values in Eq. (31) for a range of values of /3e, 
/3e w and L. 


VI. DROPLET PROFILES OBTAINED USING 
THE BINDING POTENTIAL 


Using the binding potential calculated from the DFT, 
a drop thickness profile can be calculated from the IH 
model, by minimising the free energy in Eqs. ([ 2 ]) or (|3). 
This drop profile, despite being the result of a mesoscale 
calculation, contains information about the nature of the 
microscopic interactions between particles in the system 
via the binding potential. Of course, a drop profile can 
also be calculated directly using DFT; the result of such 
a calculation is shown in Fig. [ 6 ] However, it is computa¬ 
tionally easier to treat larger systems using the IH model. 
Also, non-equilibrium situations are much more easily 
modelled via Eq. Q than with a dy nami cal DFT model 
that includes all the hydrodynamics!^^ However, since 
the two approaches are both based upon the same micro¬ 
scopic interactions, the resulting drop profiles from each 
method should be the same at the mesoscopic scale. This 
section shows how the two approaches compare. We find 



hjcT 

FIG. 10: The binding potential for a fluid with /3e = 0.9 
and /3e w = 0.55 with truncated interaction ranges of 5 <j and 
lOOcr, for the fluid-fluid and wall-fluid interactions, respec¬ 
tively. The data calculated from the DFT model is shown 
with the fitting function Eq. ( |31| ) (parameter values are given 
in Appendix A) and the sharp-kink binding potential for the 
same parameters. The inset shows the same data on a log-log 
scale. 


that overall, the drop shape profiles from the DFT and 
from the IH model are in good agreement, a result which 
a-priori is not obvious, if bearing in mind the degree of 
coarse-graining in going from a density profile to a film 
height profile. The good agreement between the two is 
thus more than a mere consistency test. 

Performing the minimisation described in Sec. |III| for 
the lattice-gas DFT, constrained via Eq. ([27]) , on a 2D 
domain yields density profiles such as those displayed in 
Fig.| 6 ] These 2D profiles correspond to the density pro¬ 
file of a cross section through a liquid ridge on a surface 
with centre of mass along the line x = 0. The loca¬ 
tion of the liquid-gas interface can be calculated using 
Eq. |7j. The alternative procedure, which yields a very 
similar result, is to just plot the contour where the den¬ 
sity pi = ( p g + pi)/ 2 = 0.5. Note that this treats the 
(strictly) discrete density profile as a continuous func¬ 
tion. In Fig. El are displayed drop profiles obtained from 
DFT and using Eq. Q, compared with results obtained 
from minimising Eq. (pj) together with the binding poten¬ 
tial obtained from DFT for various values of /3e w . These 
results show that drop profiles obtained via the two meth¬ 
ods coincide well with each other, as was also observed 
in Ref. S3 In Fig. [12] are comparisons of a droplet pro¬ 
file found via the DFT route with droplets calculated 
from a minimisation of both the full curvature free en¬ 
ergy Eq. <© and the long wavelength approximation free 
energy, Eq. ©. This shows that the full curvature free 
energy ([ 2 ]) gives closer agreement with the DFT results 
than the approximation in Eq. <©> as one would expect. 
Note too that when the contact angle is very small, then 
the drop shape is very sensitive to variations in the pa- 
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FIG. 11: Comparison of drop profiles calculated from DFT 
(points) and from the IH model, Eq. § (dashed lines). We 
compare drops of equal maximum height for a range of values 
of the wall attraction strength t w . The two methods coincide 
well across the range of contact angles. These are for the 
case when the particle interactions are truncated to a range 
of L = 5 with a strength of /3e = 0.9. The parameters used in 
the fitting functions for g(h), used to obtain the IH results, 
can be found in Appendix A. 
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FIG. 12: Drop profiles calculated from the DFT density pro¬ 
files together with those from the IH model using both the full 
curvature free energy and the long wavelength approxi¬ 
mation, Eq. These results are for /i c — 0.9 and /ie,,, — 0.5 
(top figure) and /3e w = 0.7 (bottom figure). In the inset we 
show a magnification of the contact line region. Note the ‘pre¬ 
cursor’ film height a indicating that this is sub-monolayer. 


rameters. 

Determining the contact angle from drop profiles ob¬ 
tained using the IH model is not straight-forward because 
the profiles have a precursor film. The DFT calculations 
show (see e.g. Figs.[l][5]and|6]) positive adsorption on the 
surface outside of a droplet. However, the fluid density in 
the first (i = 1) layer of lattice sites is substantially less 
than one - see also the inset in Fig. [12] This indicates 
that the density on the surface outside of the drops is not 


even a complete monolayer and so calling it a “precursor 
film” is arguably misleading. 

Returning to the issue of calculating the contact angle 
in a IH pre-cursor film model: here, this is calculated 
from the curvature at the top of the drop, at x = 0. Liq¬ 
uid droplets of a suitable size have the shape of a spherical 
cap. This is the case when the volume of the drop is small 
enough that gravity does not play a significant role but 
also large enough that the shape at the maximum is not 
distorted by the binding potential. It is droplets of this 
size that are studied here. In these cases where the drops 
have this spherical cap shape, the contact angle for the 
drop is calculated by fitting a circle to the highest point 
of the drop profile. The contact angle that this circle 
makes with the substrate is taken to be the contact an¬ 
gle of the drop. This procedure is straight-forward: The 
curvature at a stationary point, i.e. the maximum point 
of the liquid drop, is simply the second derivative of the 
height profile at that point and so the radius of curvature 
is given as 


•Emax j 

where x max is the point where the height of the droplet h 
is at it’s maximum value. Defining a to be the distance 
from the point that the circle meets the substrate to the 
centre of the base of the drop, then the contact angle is 
given by 

cos (£) - (34) 

Comparisons made in this way can only be used to 
relate the two theories for systems where a liquid drop 
exists, i.e. for a partially wetting system. The focus here 
is also on drops with an acute contact angle as these pro¬ 
files can be calculated by minimising Eq. ([3] , but droplets 
with obtuse contact angles cannot. For non-equilibrium 
problems one makes the assumption of small contact an¬ 
gles (long-wavelength) and utilises the thin film equation, 
Eq. Q, with the free energy functional given by Eq. ([3]. 
We should also mention that making the long-wave ap¬ 
proximation, i.e. going from Eq. © to Eq. (|3j) , also re¬ 
sults in the drop shape away from the surface ceasing 
to be the arc of a circle and instead being a parabola.^ 
Thus, for the long-wave theory §, we fit the drops with 
a parabola, in order to make a fair comparison. Fig. |13| 
shows the comparison of contact angles calculated for 
different values of (3e w , from the DFT model and from 
the IH model, for both the full curvature case and the 
long-wavelength approximation. 

An alternate way to obtain the contact angle is by cal¬ 
culating the largest gradient of the drop profile. The in¬ 
verse tangent of this gives the contact angle. These ‘max¬ 
imum gradient’ contact angles are also shown in Fig. [13] 
This figure shows that results based on the full curva¬ 
ture free energy agree well with those calculated from 
the DFT model. The IH model in the long-wave ap¬ 
proximation only agrees with the DFT results for very 
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FIG. 13: The contact angle for a range of different values 
of /3e w , for /3e = 0.9, calculated from both the DFT and IH 
models. Contact angles are calculated for both the full curva¬ 
ture and long-wave approximation versions of the IH model, 
using both the circle fitting and maximum gradient methods. 
The circle fitting method does not well represent the actual 
contact angle of the droplet from the long-wave IH model. 


small contact angles, as expected. For the full curvature 
model, the contact angle found by fitting a circle to the 
profile matches the DFT results better than those from 
the maximum gradient method. The droplet shape cal¬ 
culated under the long-wavelength approximation is no 
longer that of a spherical cap, it is instead parabolic, 
and so a similar procedure to the circle fitting method is 
employed where a parabola is fitted to the droplet pro¬ 
file. Results for the contact angle obtained via these fits 
are shown in Fig. [13] Clearly, these results depend heav¬ 
ily on the method used to find the contact angle of the 
droplets and although the circle/parabola fitting method 
gives better results it is only applicable to specific drop 
shapes - see Ref. [22] for a further discussion on extracting 
a contact angle from a drop profile. Note also that the 
contact angles in Fig. |T3] calculated via DFT are the only 
results to extend beyond 90°; droplet profiles can not be 
found in this range using the IH model. 

The calculated contact angles are weakly dependent on 
the size of the drop. All contact angles extracted from the 
IH model are calculated from droplets that have a height 
of 35cr. Calculating contact angles from droplets of a spe¬ 
cific height, rather than a specific volume, seems to give 
more consistent results. The range at which particle in¬ 
teractions are truncated does not significantly affect how 
well the results from the two models (DFT an IH) agree 
with each other. Over the range of parameters studied, 
up to a truncation range of L = 40, the discrepancy in 
the contact angle between the two models is typically a 
few percent. 


VII. CONCLUSIONS 


In this paper we have developed a microscopic DFT 
based method for calculating the binding potential g(T) 
at various different state points and interaction poten¬ 
tials of various different strengths and ranges between the 
surface and the liquid. These were subsequently used as 
input to a mesoscopic interface free energy model, which 
was then used to calculate the height profile of liquid 
drops on surfaces. The liquid height profiles are very 
similar in shape to the profiles obtained directly from 
the DFT, indicating that the coarse-graining procedure 
used to obtain g(T) is valid. 

The binding potential is calculated as a function of the 
adsorption T, which can be related to the film height via 
Eq. ( fl3|). It is calculated using the method developed 
in Ref. 1471 for studying nucleation of the liquid from the 
gas phase. This method constrains the fluid density pro¬ 
file to have a specified value of T, which is equivalent to 
imposing an additional fictitious external potential that 
stabilises the adsorbed him with the given adsorption. 
Note that this fictitious potential is not known a-priori 
and is calculated on-the-fly as part of the minimisation 
to obtain the fluid density profile. 

The method presented for calculating the binding po¬ 
tential is general and it should be possible to use if for 
any DFT. Here it has been implemented using a simple 
DFT for a lattice-gas, that in Eq. (23). A better descrip¬ 
tion of the free energy of a lattice-gas is possible, we refer 
the interested reader to e.g. Ref. E3 for a more sophisti¬ 
cated and accurate approximation for the lattice-gas free 
energy. The lattice-gas approximation was used in or¬ 
der to more easily test the results for a range of different 
drop sizes. Note that to calculate the 2D density pro¬ 
files for drops of the size in Fig. [6] (or even larger) with 
a more sophisticated continuum DFT can be computa¬ 
tionally time consuming - e.g. using the commonly used 
approximation F[p{ r)] = J : fmt[p(*)\+ F a tt[p{*)\, where 
Ffmt[p( r)] is the fundamental-measure theory approxi¬ 
mation of the free energy of a fluid of hard-spheres and 
T a tt\p{ r)] is a simple mean-field approximation for the 
contribution to the free energy due to the attractive in¬ 
teractions between the fluid particles ISUUtHI However, one 
of the draw-backs of the lattice-gas model is that the dis¬ 
cretization onto a lattice leads to unrealistic small ampli¬ 
tude oscillations in g(T ), particularly for larger values of 
/3e (i.e. for lower temperatures) - see e.g. Fig. [7j A more 
realistic continuum DFT model will not exhibit such os¬ 
cillations that do not decay in amplitude with distance 
from the wall. Note, however, that oscillations may be 
present in the true binding potential. These, if present, 
stem from the packing of the particles, rather than from 
the discretization onto a lattice and are to be expected 
whenever the fluid exhibits layering transitions. 

Another drawback of the lattice-gas is that the map¬ 
ping of the continuum fluid onto the lattice constitutes a 
significant approximation. This prevents a quantitative 
comparison between our results with existing simulation 
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results for continuum models, such as those in Refs. HU¬ 
ES] However, since we do find good qualitative agree¬ 
ment, we are currently applying the method using con¬ 
tinuum (FMT) DFT. Results from this follow-on study 
will be published elsewhere, including comparison with 
simulation results. 


It should also be emphasised that it remains straight¬ 
forward to calculate droplet profiles using DFT rather 
than using the mesoscale model over the full range of 
contact angles, i.e. including when the contact angle > 
90°. However, the IH model is particularly advantageous 
when considering the non-equilibrium situation, which 
can be described using Eq. Q, albeit this equation is 
derived under the assumption of small contact angles and 
only accounts for convective transport with no slip at the 
substrate. It is shown in Fig.[l3]how the accuracy rapidly 
breaks down for larger contact angles. Non-equilibrium 
phenomena can also be studied using dynamical DFT, 
see e.g. Refs. fSBHSTI for more on this approach. 

One of the striking results of the present work is the 
observation of the degree to which the binding potential 
changes when the range at which particle interactions 
are truncated is changed - see in particular Figs. [ 8 ] and 
[9j Crucially, the value of g(r) at the minimum that can 
be present at low values of the adsorption depends very 
sensitively on the truncation range. This minimum can 
change from being just a local minimum (a metastable 
equilibrium point) to become the global energetic mini¬ 
mum which represents a shift in the phase behaviour from 
wetting to non-wetting. In fact, for very short truncation 
ranges, the minimum can be completely removed. The 
conclusion from this part of our study is that to deter¬ 
mine accurately the location of a wetting transition in 
theory or simulations one must ensure that the trunca¬ 
tion is not so severe as to induce such errors. The tails 
of the potentials matter! 

This is important in the context of (coarse-grained) 
Molecular Dynamics simulations where dispersion forces 
are often cut off at distances of a few particle diame¬ 
ters. Ref. [ 22 ., for instance, employs a cut-off length for 
the Lenard Jones interactions of two times the equilib¬ 
rium bead distance and finds, in consequence, that the 
extracted binding potential is well fitted by a sum of ex¬ 
ponentials as in Eq. similarly to our results when 

only short-range interactions are considered. 

Note also that for weakly attracting (solvophobic) 
walls the minimum in g(T) is at very small, or poten¬ 
tially even negative, values of T. In this regime describ¬ 
ing the sub-monolayer adsorption via a film-thickness h 
is a somewhat misleading concept. In this regime, the 
minimum in the binding potential is very asymmetric, 
rising very sharply on the small T side of the minimum, 
but rising far less steeply as T is increased from the value 
at the minimum. To incorporate this behaviour in the 
fit-function for g(T), we had to use the exponential in 


the function P(T) in Eq. (31). This also means that the 
minimum in g(r) is not we 1 -approximated by a quadratic 
function when the minimum is at a small value of T and 


so, of course, capillary-wave theory also does not apply 
in this regime, since capillary-wave theory assumes Gaus¬ 
sian fluctuations in a harmonic potential. Also, when the 
adsorption at the surface is sub-monolayer, then a non¬ 
equilibrium situation can no longer be described by Eq. 
Q; it can instead be described via a gradient dynamics 
model with a diffusive dynamics. 

Finally, it should also be pointed out that although in 
the present work we have developed a method for calcu¬ 
lating a better approximation for the binding potential, 
our approach is still a mean-field theory and therefore 
does not include all fluctuation effects, because our the¬ 
ory is based on Eq. §■ Some effects of fluctuations at 
wetting transitions can be included by replacing 7 i g in 
Eq. ([ 2 ]) with a function 7 (h). However, as Parry and co¬ 
workers have showed, the effective interface Hamiltonian 
is in fact non-localP^HMl Non-locality and fluctuation ef¬ 
fects are particularly important at wetting transitions. 


Acknowledgements 

A.P.H. acknowledges support through a Loughborough 
University Graduate School Studentship. All authors 
thank the Center of Nonlinear Science (CeNoS) of the 
University of Munster for recent support of the authors 
collaboration. 


Appendix A: Fit function parameters 


The following table contains the parameters to the fit¬ 
ted binding potentials that have been used throughout 
this paper and also for a selection of other state points. 
For ease of reference, the fitting function is repeated here: 


sin = 


(Al) 


with 

P(r) = a 0 r 2 e" air + a 2 r 2 + a 3 r 3 + a 4 P 4 + a 5 r 5 . 


Appendix B: Particle interaction potentials 


The net interaction potential between a fluid particle 
and a planar wall made of particles interacting with the 
fluid particle via a Lennard-Jones-like potential, with the 
same form as in Eq. (21) is 


0 00 00 

Vi = e wf Y. E E (Bl) 

i '=—00 j '=—00 k '=—00 


where we assume the surface of the wall is located in 
the plane * = 0. The parameter e w f characterises the 
strength of the interaction between a single wall particle 
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TABLE I: This table lists the parameters used in the fitting function to generate the results shown in the figures within the 
main text. The parameters are identified by the figure in which the binding potential is used, the title of that curve in the 
figure and some additional identification where applicable. Values are rounded to three decimal places and the value of a\ is 
enforced. 


and a fluid particle. The triple sum above spans all of 
the particles in the wall; i.e. the wall is modelled as being 
discretised on a lattice just as the fluid is and that a wall 
particle is present on every lattice site where i < 0. The 
sums in Eq. (B11 can be simplified to obtain 


Vi = 


-e w /i 3 


for i > 1 
for i < 1 


(B2) 


The parameter e is the strength of a single Lennard-Jones 
pair interaction between two fluid particles and eij is the 
weighted interaction between two particles taking into 
account all of the interactions in the invariant k dimen¬ 


sion. In the right hand side of Eq. (B3), i', j’, and k! are 


the distances between a pair of lattice sites in the i, j 
and k directions respectively. 


where the parameter e w determines the net strength of 
attraction to the wall. 

The (2D) effective fluid-fluid particle interactions are 
governed in a similar manner by the potential 


e U 


e (i’ 2 + j' 2 )- 3 + 2 Y, (*' 2 + f + k ' 2 )~ 3 

k ’=1 


(B3) 


Due to the fact that calculating long ranged particle 
interactions can be time consuming, making computer 
calculations rather slow, it is often necessary to truncate 
the interactions to some interaction range La. When 
particle interactions are truncated at a range of La the 
external potential Eq. (B11 becomes 


Vi = e 


o 


E E 


i'——L j'——L k'——L 


(( i-i'f+f 2 


/ 2 \ —3 


k’ 2 ) 


for ((a - i’f + j' 2 + k’ 2 ) < L 2 
otherwise 


(B4) 
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and the fluid-fluid particle interaction weights (B31 are then given as 



f 

k'<y/L 2 -i' 2 -j' 2 


e 

(i' 2 +j' 2 )~ 3 + 2 £ {i' 2 + j' 2 + k' 2 )~ 3 

n 

,o 

k '=1 


for |i — j| < L 
otherwise. 


(B5) 
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Note that the contribution to the lattice ga s fr ee energy 
density /3f = plnp — (1 — p) ln(l — p) in Eq. (231 contains 
both the ideal gas /3 fid = plnp — p and an excess contri¬ 
bution 0f ex = p + (1 - p) ln(l - p); i.e. f = fid + fex- The 
origin of this form arises from the hole-particle m <->• 1 — m 
symmetry in the lattice model (see e.g. Ref. 14511 . The sec¬ 
ond term, f ex , stems from the particles having a hard-core 
interaction, preventing multiple occupancy of a cell. It is 
instructive to compare this with the free energy of a bi¬ 
nary mixture. The ideal-gas contribution to the free en¬ 
ergy of a mixture is fif = pi hi pi — pi + p 2 \np 2 — p 2 , 
where pi is the density of species 1 and p 2 the density 
of species 2. Re-writing this in terms of the total den¬ 
sity p = pi + p 2 and concentration x = pi/p, we obtain 
/3f = p In p — p + p[x In x + (1 — x) ln(l — x)\. The final term 
is the ideal-gas contribution to the entropy of mixing and 
it is interesting to note the similarity in form to the lattice- 
gas free energy. For an incompressible binary mixture p is a 
constant and so the first two terms are often neglected. As 
a consequence of this, the free energy of an incompressible 
binary fluid can have a very similar mathematical form to 
the lattice-gas free energy considered here. 



